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Abstract . A pressurized cylindrical shell containing a 
longitudinal crack is considered. It is assumed that the 
material has a special orthotropy, namely that the shear modu- 
lus of the sheet may be evaluated from the measured Young's 
moduli and the Poisson's ratios rather than being an independ- 
ent material constant. Two examples, one for a mildly ortho- 
tropic (titanium) and the other for a strongly orthotropic 
(graphite) material approximately satisfying the condition of 
special orthotropy are given. The results show that the stress 
intensity factors are rather strongly dependent on the degree 
of orthotropy. 


SOME REMARKS ON THE FORMULATION OF THE PROBLEM 


A detailed treatment of the linear bending theory of aniso- 
tropic shallow shells may be found in [1-3]. Assuming that 
through a proper superposition the problem has been reduced to 
a perturbation problem with the crack surface tractions as the 
only external loads, in an eight order theory the differential 
equations for an orthotropic shallow cylindrical shell may be 
expressed as 


D l v l w ( x i»x 2 ) R 9 2 F ( x i > x 2 ) 

o X i 

a h E p s g 2 

V 2 F ( x i » x 2 ) pi 2 ^ ( x i » x 2 

9 x -j 


= 0 , 


« o , 


(1 .a,b) 
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where a,R,h are the dimensions of the crack and the shell (see 
the insert in Figure 1), 

x-j = X/a , x 2 = Y/a , (2) 


F is a stress function and w is the displacement component nor- 

4 4 

mal to the surface. The operators V-j and v 2 are defined by 


^4 ^1 2 

+ 2[v 0 + 2(1-v ] v 2 ) - 4 1 ] 


3x 


3 x 


+ 2 ( 


1 


3 x -j 3 x 2 
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1 3 X 


2G 


12 


- v 2 ) 


3x-j 3x 2 


2 + E 


1 3x 


4 * 


(3. a ,b) 


The notation for the orthotropic elastic constants are defined 
by the following stress strain relations: 

e ll = Fj" l°ll ' V 1 a 22 ^ ’ e 22 = ( ct 22 " v 2 a ll ^ ’ 

£ 1 2 = 2G^ a 12 ’ = if ’ D 1 = E-j h 3 / [12(1 -v-j v,,)] , 

(4 . a-e) 

where x-j and x 2 are the principal directions of orthotropy and 
are taken respectively along the axial and circumferential 
directions. The stress and moment resultants are related to 
F and w through the usual expressions (see [3] or [4]). 


In solving the problem, for example by expressing F and w 
in terms of appropriate Fourier integrals, (1) may be reduced 
to a system of two fourth order ordinary differential equations. 
The characteristic function of this system will then be an 8th 
degree polynomial the coefficients of which will be functions 
of the transform variable. For the problem to be analytically 
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tractable it is essential that the roots of the characteristic 
equation be obtainable in closed form. Apparently, for the 
anisotropic shells in general and for the orthotropic shells 
in particular this is not possible. In order to express the 
roots in closed form the operators V-j^ and must be properly 
factorized. From (3) it is clear that these operators can 
indeed be factorized and may be expressed in the following form 

= (A + ^A > 2 = v • < 5 > 

o x i o X ^ 


provided the elastic constants satisfy the following conditions: 

G 


[v 2 + 2(1-v 1 v 2 ) = 1 

' v 2^ ’ /E 1 /E 2 = 1 ' 


( 6 . a , b) 


By direct substitution it may easily be shown that the conditions 
(6) are satisfied if 


G 


12 


( E l E 2^ 

2 [1 + (v 1 v 2 )' 5 ] 


(7) 


Considering also the relation (4.d), this means that the sheet 
material has only three independent constants. Such a material 
is said to be specially orthotropic. The results given in this 
note will then be valid only for those materials in which the 
measured value of G ]2 and that calculated from (7) in terms of 
measured and v - , (i=l,2) are in reasonably good agreement. 


Changing the variables once more as 
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( 8 ) 


X = X 


(5) becomes 


(e 1 /e 2 )^ x 2 = y 
2 2 


V , 4 ■ v 4 = (^-y. + ^-g-) = v‘ t . 

8x 3y 


(9) 


With (9), (1) becomes identical to the differential equations 

3 p 

for isotropic shells in which D = Eh /[I 2(1 -v )] and E are 
replaced by D-j and respectively. The problem may then be 
solved by following, for example, the same procedure as that 
described in [4] or [5] and will not be repeated here. 


STRESS STATE AROUND THE CRACK TIP AND STRESS INTENSITY 

FACTORS 

Following very closely the technique described in [4], for 
a pressurized shell the membrane and bending components of the 
stresses in the close neighborhood of the crack tip may be 
obtained as 

a vv V,e) = A m — jj; ^ (5cos|- - cosrjr) + O(r^) , 
y y /2 r 

a xx m (r, e ) = A m — ^ Sp (3cos| + cos^-) + O(r^) , 

o, *(r.8) = A — x (sin| - sin^-) + Ofr* 4 ) , (lO.a-c) 
xy m ^ Hi i 


k 5 v v 

a yy (r,0sZ) = A b ^ 2F [(3v o ' 8v c + _ ^P' )c ° s f 


V 0 V C „ 50 


2 ^ cos^] + 0(r*) 
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b/ - ~ - fl —2- — r ( ft + — £ - 8v c 2 + 3v 0 v )cos|- 

2 C J 


xx (r.e.z) - A b [(8 + -2 


4c 


( 1 -v 9 c 2 ) cos^|-] + 0(r" 2 ) 


b, , _ . k p z a2h3(3 i 2 r v c . . 50 

°xy (r>e>z) = A b h ~T2D^ — [ ^2 sln T 


- ( 8 + -&) s 1 ni] + 0 ( r" 5 ) , 


(11 .a-c) 


where (r,9) are the polar coordinates at the crack tip, z is the 
coordinate perpendicular to and measured from the neutral surface 
of the she! 1 , and 


p R _ 

■y/a , c = ( E ] / E 2 ) 


v 0 = c - V, 


_ 2 , h|3G 12 1 
v c - c ' ^ V 1 3D„ ^ ’ 


E ^ h 

1 2 ( 1 - v-j Vg ) • 


( 12 ) 


Here p Q is the internal pressure. Thus, kp is the corresponding 
stress intensity factor in the flat plate, and A m and A b repre- 
sent stress intensity factor ratios defined by 


A. = k*/k 


m 


m / “p 5 


A b ■ k b /k P 


(13. a, b) 


where k s and k? are the membrane and bending components of the 
mb 

stress intensity factor in the shell. 


NUMERICAL EXAMPLES 

The elastic constants of the orthotropic shells which are 
considered as examples are shown in Table 1. The Table also 
shows the "average shear modulus" calculated from 
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(14) 


( Ei e 2 )^ 


av 


2 [ 1 + (v-j v 2 ) 2 ] 


Table 1. Elastic Constants of the Materials 



Titanium 

Graphite 

E] (psi) 

1 ,507xl0 7 

1 . 5x1 0 6 

E 2 (psi) 

2 . 08x1 0 7 

40x1 0 6 

v i 

0.1966 

0.0075 

v 2 

0.2714 

0.2000 

G 12 

6.78x10® 

4.0x10® 

G 

7. 15x10® 

3.73x10® 

av 




Table 1 indicates that for the two materials under consid- 
eration G 10 and G are sufficiently close to justify the 
12 av 

assumption of special orthotropy. The calculated stress inten- 
sity factor ratios based on this assumption are shown in Tables 
2 and 3 and Figures 1-4. In the examples E-j is taken in the 
axial direction. The tables also show the related (orthotropic) 
shell parameter X Q calculated from 

X 0 = [12(1-v 1 v 2 )E 2 /E 1 ] % a/tRh)^ . (15) 

Since X Q is dependent on the elastic constants, in comparing 
the orthotropic results for different material orientation and 
the results for the isotropic material, a/(hR)^ is used as the 
independent variable. 
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Table 2. 


The Stress Intensity Factor Ratios 
for a Pressurized Titanium Cylinder 


(E 1 /Eg) = 0.724 


X o 

a/(Rh) % 

A 

m 

A b 

1 

0.5025 

1.1962 

0.1453 


1 .0050 

1 .6053 

0.2645 


1.5075 

2.0795 

0.3091 


2.0100 

2.5638 

0.2853 


2.5126 

3.0419 

0.1983 


3.0151 

3.5094 

0.0590 


3.5176 

3.9648 

-0.1 255 


4.0201 

4.4081 

-0.3493 

( E-j /Eg) = 1.380 

1 

0.5903 

1.1962 

0.2010 

mm 

1 .1807 

1.6053 

0.3643 

B 

1.7711 

2.0794 

0.4265 


2.3615 

2.5638 

0.3927 


2.9518 

3.0420 

0.2734 

B 

3.5422 

3.5094 

0.0808 


4.1326 

3.9649 

-0.1743 

H 

4.7229 

4.4082 

-0.4835 
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Table 3. The Stress Intensity Factor Ratios 
for a Pressurized Graphite Cylinder 


(E-j /E 2 ) = 0 . 0375 



a/ ( Rh) ^ 

m 

A b 

1 

0.2365 

1.1928 

0.0217 


0.4730 

1 . 5906 

0.0464 


0.7095 

2.0493 

0.0645 

n 

0.9460 

2.5179 

0.0726 


1.1826 

2.9818 

0.0697 

n 

1 .4191 

3.4373 

0.0564 

n 

1 . 6556 

3.8833 

0.0338 

n 

1 .8922 

4.3196 

0.0027 

( E -j / E 2 ) = 26.66 

1 

1 .2214 

1.1928 

0.5787 

mm 

2.4428 

1 . 5906 

1 .2393 


3.6642 

2.0493 

1.7220 

D 

4.8856 

2.5179 

1 . 9368 

WM 

6.1070 

2.9818 

1.8593 

El 

7,3284 

3.4373 

1.5049 

n 

8.5498 

3,8833 

0.9007 

Hi 

9.7712 

4.3196 

0.0726 
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For the purpose of comparison, Figures 1-4 also show the 
results for an isotropic shell (with a Poisson's ratio of 1/3) 
which are given in [5], It is seen that in the specially ortho- 
tropic shells the stress intensity factors are strongly dependent 
on the modulus ratio E i / E ^ > and A m and A^ generally increase with 
decreasing E i / E 2 • However, this does not necessarily mean a 
reduction in the resistance of the shell to cleavage as E -j / E 2 
decreases (i.e., as the shell becomes stiffer in circumferential 
direction). For this, one also has to consider the fracture 
strength of the shell in the plane parallel to E-j as a function 
of E-j/Eg. Intuitively, it is expected that this strength too 
would increase as E 1 / E 2 decreases. 
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Figure 2. Bending component of the stress intensity 
factor ratio for a pressurized Titanium cylinder. 



Figure 3. Membrane component of the stress intensity 
factor ratio for a pressurized Graphite cylinder. 



